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Abstract. The origin of the extragalactic gamma-ray background (EGRB) radiation has been 
a mystery in astrophysics for a long time. Recently the Fermi gamma-ray satellite (Fermi) has 
revealed that ~22% of the unresolved EGRB would be explained by blazars, which are one 
population of beamed active galactic nuclei (AGNs). The remaining ~78% of the unresolved 
EGRB is still unknown. We estimate the contribution of gamma-ray loud radio galaxies, which 
are misaligned radio loud AGNs recently detected by Fermi, to EGRB using the radio luminosity 
function of radio-loud AGNs with the correlation between the radio and gamma-ray luminosities. 
We find that ~25% of the unresolved EGRB will be explained by gamma-ray loud radio galaxy 
population. We also discuss further about the origin of EGRB by comparing the Fermi EGRB 
data with our studies on various AGN populations' contribution to EGRB, which are radio 
quiet AGNs, blazars, and gamma-ray loud radio galaxies. 



1. Introduction 

The origin of the extragalactic gamma-ray background (EGRB) radiation has been argued for 
a long time in astrophysics. Although several sources have been suggested to explain the MeV 
and GeV background, it is still discussing about the origin of the MeV background due to the 
difficulties of MeV gamma-ray measurements. In the case of the GeV background, since blazars 
which are a type of active galactic nuclei (AGNs) are dominant extragalactic gamma-ray sources 
[TJ |2], it is expected that unresolved population of blazars would explain the GeV EGRB (see 
e.g. [3]). Recently, [4] showed that unresolved blazars can explain ~ 22% of the EGRB above 0.1 
GeV by analyzing the 11-months Fermi AGN catalog. Other gamma-ray emitting extragalactic 
sources have also been discussed as the origin of the GeV EGRB (see [5] for details). 

Fermi has recently detected GeV gamma-ray emissions from 11 misaligned AGNs (i.e. radio 
galaxies), which are one type of AGNs with the direction of a relativistic jet not coinciding with 
our line of sight [6]. Their average photon index at 0.1-10 GeV is ~ 2.4 which is same as that 
of GeV EGRB [7] and blazars [3]. Although they are fainter than blazars, the expected number 
in the entire sky is much more than blazars. Then, it is naturally expected that they will make 
a significant contribution to EGRB. 

In this proceeding paper, we discuss the contribution of gamma-ray loud radio galaxies 
(not blazars) to EGRB based on [5]. Furthermore, we also discuss the EGRB contribution 
from various AGN populations such as radio-quiet AGNs and blazars based on our previous 
studies \19\ [3]. Throughout this paper, we adopt the standard cosmological parameters of 
(h,n M ,n A ) = (0.7,0.3,0.7). 



2. Gamma-ray Luminosity Function of Gamma-ray Loud Radio Galaxies 

To estimate the EGRB contribution of gamma-ray loud radio galaxies, we need to construct a 
gamma-ray luminosity function (GLF). However, because of a small sample size, it is difficult 
to construct a GLF using the current gamma-ray data only. Here, the radio luminosity function 
(RLF) of radio galaxies is extensively studied by previous works (see e.g. [101 [11]). If there is a 
correlation between the radio and gamma-ray luminosities, we are able to convert the RLF to 
the GLF with that correlation. 

Figure [T] shows the 5 GHz and 0.1-10 GeV luminosity relation of the Fermi gamma-ray loud 
radio galaxies. Square and triangle data represents FRI and FRII radio galaxies, respectively. 
The solid line shows the fitting line to all the data. The function is given by 

log 10 (L 7 ) = (-3.90 ± 0.61) + (1.16 ± 0.02) log 10 (L 5G H z ), (1) 

where errors show 1-a uncertainties. In the case of blazars, the slope of the correlation between 
L 7 (> lOOMeV), luminosity above 100 MeV, and radio luminosity at 20 GHz is 1.07 ± 0.05 [8]. 
The correlation slopes of gamma-ray loud radio galaxies are similar to that of blazars. This 
may indicate that emission mechanism is similar in gamma-ray loud radio galaxies and blazars. 
The partial correlation coefficient [9] for these data set becomes 0.866 with chance probability 
1.65 xlO" 6 . Therefore, we conclude that there is a correlation between the radio and gamma-ray 
luminosities of gamma-ray loud radio galaxies. 

We derive the GLF of gamma-ray loud radio galaxies, p 7 (L 7 , z), using the correlation between 
the radio and gamma-ray luminosities as Equation [1] and the RLF of radio galaxies, p r (L r ,z), 
with radio luminosity, L r . The GLF is given as 

p 7 (L 7 , z) = K^-p r (L r , z), (2) 

where k is a normalization factor. We use the RLF by |llj . 

The normalization factor n, which corresponds to the fraction of gamma-ray loud radio 
galaxies against all radio galaxies, is determined by normalizing our GLF to the source count 
distribution of the Fermi radio galaxies, which is sometimes called logN-logS plot or cumulative 
flux distribution. Source count distribution is calculated by 

Ar(>F 7 ) = 47r/ dz—— dL^p 7 (L 7 ,z), (3) 

Jo dilaz jl 7 ( 2 ,f 7 ) 

where L 7 (z,i ? 7 ) is the gamma-ray luminosity of a blazar at redshift z whose photon flux at 
>100 MeV is F 1 . Hereinafter, we assume z max — 5 and L 7j max — 10 ^erg/s in this study. These 
assumptions hardly affect the results in this study. 

Since the completeness of the Fermi sky survey depends on the photon flux and photon 
index of a source, we need to take into account this effect (so called the detection efficiency) to 
compare GLF with the cumulative source count distribution of gamma-ray loud radio galaxies. 
The detection efficiency of Fermi is shown in Figure 7 of [3]. 

Figure [2] shows the source count distribution of gamma-ray loud radio galaxies. The data 
is after the conversion of the detection efficiency. Solid red curve shows in the case of k = 1 
which corresponds to the case that all radio galaxies emit gamma-rays. Dashed blue curve 
corresponds to the GLF fitted to the Fermi data with k = 0.081 ± 0.011. ~ 1000 gamma-ray 
loud radio galaxies are expected with 100% complete entire sky survey above the flux threshold 
F 7 (> lOOMeV) = 1.0 x 10~ 9 photons cm" 2 s" 1 above 100 MeV. We note that the current 
detection efficiency of Fermi at F 7 (> lOOMeV) = 1.0 x 10~ 9 photons cm -2 s _1 is ~ 10~ 3 . 
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Figure 1. Gamma-ray luminosity at 0.1-10 GeV 
versus radio luminosity at 5 GHz. The square and 
triangle data represents FRI and FRII galaxies, 
respectively. The solid line is the fit to all sources. 
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Figure 2. The source count distribution of 
gamma-ray loud radio galaxies in the entire sky. 
Solid and dashed curve corresponds to the all 
radio galaxies and gamma-ray loud radio galaxies, 
respectively. The Fermi data are also shown. 



3. EGRB from Gamma-ray Loud Radio Galaxies 

We calculate the EGRB spectrum by integrating our GLF in the redshift and luminosity space, 
using the a simple broken power law SED model (see [5] for details). The EGRB spectrum is 
calculated as 



rL-y, max dL[Lj, (1 + z)e] 

I aL 7 /9 7 (L 7 , z) '— exp[— r 7i7 (e, z)\ x{1.0— cj(F 7 [L 7 , z\) }, 

-^7,min 

where t is the cosmic time and dt/dz can be calculated by the Friedmann equation in the standard 
cosmology. The minimum gamma-ray luminosity is set at £ 7) min = 10 39 erg/s. w(F 7 [L 7 ,z]) is 
the detection efficiency of Fermi at the photon flux F 7 (L 7 , z) . 

Here, high energy 7-rays propagating the universe are absorbed by the interaction with the 
extragalactic background light (EBL), also called as cosmic optical and infrared background 
[T2J [13] . r 7 , 7 (e, z) is the optical depth of this background radiation. In this study, we adopt the 
EBL model by [13]. We also include the contribution of the cascade emission. The gamma-ray 
absorption creates electron-positron pairs. These pairs scatter the cosmic microwave background 
radiation to make the secondary emission component (so called cascade emission) to the absorbed 
primary emission |14j . 

Fig. [3] shows the unresolved vl v EGRB spectrum of gamma-ray loud radio galaxies in the 
unit of MeV 2 cm _2 s _1 MeV~ 1 sr~ 1 predicted by our GLF. The intrinsic (the spectrum without 
the EBL absorption effect), absorbed, and cascade components of the EGRB spectrum and the 
total EGRB spectrum (absorbed+cascade) are shown. As in the figure, the cascade emission 
does not contribute to the EGRB spectrum significantly. 

The expected EGRB photon flux above 100 MeV from gamma-ray loud radio galaxy 
populations is 0.26 x 10~ 5 photons cm~ 2 s _1 sr . As the unresolved Fermi EGRB flux above 
100 MeV is 1.03 x 10" 5 photons cm 2 s 1 sr 1 [7], the gamma-ray loud radio galaxies explains 
~ 25% of the unresolved EGRB flux. For the comparison, recent analysis of Fermi blazars 
showed that blazars explains ~ 22% of the unresolved EGRB [3]. Therefore, we are able to 
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Figure 3. The unresolved EGRB 
spectrum from gamma-ray loud ra- 
dio galaxies. Dashed, dotted, dot- 
dashed, and solid curves show the 
intrinsic spectrum (no absorption), 
absorbed, cascade, and total (ab- 
sorbed+cascade) EGRB spectrum, 
respectively. The observed data of 
HEAO-I\jjft, Swift-BAT [E], SMM 
[T7], COMPTEL Jig, and Fermi- 
lO" 2 io -1 io° io 1 io 2 io 3 io 4 io 5 io 6 LAT [7] are also shown by the sym- 
Photon Energy [MeV] bols indicated in the figure. 

explain ~ 47% of EGRB by radio-loud AGN populations. We note that, because of small 
samples and the uncertainty in the correlation between radio and gamma-ray luminosities, the 
fraction of gamma-ray loud radio galaxies in unresolved EGRB varies from ~ 10% to ~ 63%. 

4. EGRB from Beamed and Unbeamed AGNs 

We have also conducted works on EGRB from radio-quiet AGNs [19J and blazars [3J. Raio-quiet 
AGNs would significantly contribute to the MeV background by including non-thermal electron 
population in the AGN hot corona above the accretion disk |19j . Our blazar GLF model [3] is 
constructed based on the EGRET data taking into account blazar sequence |20t |2"T] . This blazar 
GLF is consistent with various Fermi data (source counts and EGRB contribution) [22\ I23j . 
Fig. [5] shows the EGRB contribution of radio-quiet AGNs, blazars, and gamma-ray loud radio 
galaxies (i.e. beamed and unbeamed AGN populations). As in the figure, AGN populations are 
dominant in the cosmic background radiation from X-ray to gamma-ray. In this figure, we use 
the total (resolved + unresolved) EGRB data to avoid the instrument dependence, since the 
model by [3] is based on the EGRET data. 

5. Origin of the GeV EGRB 

Although AGN populations dominate EGRB, it is difficult to explain 100% of EGRB. Fig. [5] 
shows the unresolved gamma-ray loud radio galaxy EGRB spectra at each redshift bins. Because 
of EBL, the spectrum above > 30 GeV shows the absorbed signature. Here, the cosmological 
sources have their evolution peaks at z = 1 ~ 2 such as cosmic star formation history and AGN 
activity (see e.g. |24[ 125]), This means that the gamma-rays from extragalactic sources (e.g. 
galaxies and AGNs) will experience the EBL absorption. However, as shown in Fig. [SJ the Fermi 
EGRB spectrum does not clearly show such an absorbed signature. This might suggest that 
nearby gamma-ray emitting sources or sources with very hard gamma-ray spectrum would be 
the dominant population of EGRB above 10 GeV such as high-frequency-peaked BL Lacs. To 
address this issue, we should await the EGRB information above 100 GeV by future observations 
such as Fermi. CTA would also be possible to see the EGRB at much higher energy band. We 
also need to examine the EBL models at high redshift. It is expected that CTA will see blazars 
up to z ~ 1.2 at very high energy band > 30 GeV |22j . Therefore, Fermi and CTA will be a key 
to understanding the origin of EGRB. 
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Figure 4. EGRB spectrum from various AGN populations. Magenta, green, blue, and black 
curves show the EGRB contribution from radio-quiet AGN [19J, blazars [3], gamma-ray loud 
radio galaxies [5], and total AGN populations (radio-quiet AGNs, blazars, and gamma-ray loud 
radio galaxies), respectively. The observed data are same as Fig. [3l but total (resolved + 
unresolved) EGRB data are shown here. 



6. Conclusion 

In this paper, we studied the contribution of gamma-ray loud radio galaxies to the EGRB by 
constructing their GLF. First, we explored the correlation between the radio and gamma-ray 
luminosities of gamma-ray loud radio galaxies which are recently reported by Fermi. We found 
that there is a correlation L 7 oc L^q^ by a partial correlation analysis. 

Based on this correlation, we defined the GLF of gamma-ray loud radio galaxies using the 
RLF of radio galaxies. We normalized the GLF to fit to the cumulative flux distribution of 
Fermi samples. Then, we predicted the contribution of gamma-ray loud radio galaxies to the 
MeV and GeV EGRB. We found that gamma-ray loud radio galaxies will explain ~ 25% of the 
EGRB flux above 100 MeV and also make a significant contribution to the 1-30 MeV EGRB. 
Since blazars explain ~ 22% of EGRB, ~ 47% of EGRB would be explained by blazars and 
gamma-ray loud radio galaxies. 

Based on our previous studies of radio-quiet AGNs and blazars, we also showed the 
contribution of beamed and unbeamed AGN populations to EGRB. AGN populations dominate 
the cosmic X-ray and gamma-ray background in about 7 orders magnitude of photon energy. 

We also make an interpretation on the origin of the EGRB above 10 GeV from the point of 
view of the EBL absorption effect. If the EBL absorption signature is still not appeared in the 
EGRB spectrum, the origin would be nearby sources or sources with hard gamma-ray spectrum. 
We should await the EGRB data at higher energy band for this issue. 
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Figure 5. Gamma-ray loud radio 
galaxy EGRB spectra at each redshift 
bins. Solid, dashed, dotted, and 
dot-dashed curves shows the EGRB 



spectrum at z = 0.0 — 5.0, z = 0.0 — 0.8, 
z = 0.8 - 2.0, and z = 2.0 - 5.0, 
respectively. The thin and thick curve 
corresponds to the intrinsic spectrum 
(not absorbed) and the total spectrum 
(absorbed + cascade), respectively. The 
unresolved EGRB data of Fermi-LAT [7 
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